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ABSTRACT 

We derive analytical expressions for the mass flow rates in and from accretion discs, 
taking into account the Eddington limit. This allows us to connect the basic properties 
of outflows with those of the accretion flow. As an example, we derive the radial surface 
density distribution in the accretion disc of Mrk 231. 
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1 INTRODUCTION 

It is now widely accepted that active galactic nuclei (AGN) 
are powered by accretion of gas and dust via a disc on to 
the AGN's central black hole. This process releases suf- 
ficient energy to explain the obser ved AGN luminosities 
l|Lvnden- Belli Il969l : IZeTbovichl 1 1964) . At high luminosities 
(radiation) pressure may become relevant and limit the ac- 
cretion rate to the so-called Eddington limit 
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( Eddi ngton! Il92ll ). with the black h ole mass Mbh , the 
Salpeter time-scale rs = r/ ■ 4.5 • 10 8 yr jSalpeterlll964l ) and 
the accretion efficiency rj, which is ~ IP -1 for standard ac- 
cretion discs (jShakura fc Sunvae vl ll973l ). 

Thus the black hole cannot swallow all the mass that 
is delivered by the accretion disc, what happens to the 
excess mass has been a su bject of some intense debate. 
iDuschl fe Strittmatterl |201ll ). for instance, assume that this 
mass is driven from the AGN . Such outflows were i ndeed 
observed from q uasars (e.g., iMorganti et al.l (|2005l ) and 
Irlolt et all (12011 ), and, in parti cular for the context of this 
paper, iRupke fc Veilleuxl (|201ll ) from Mrk 231). 

But some aspects of these outflows remain elusive, e.g., 
what amount of mass they contain and what their power 
source is. In this paper we present analytical and numeri- 
cal investigations of mass flows in accretion discs and mass 
outflows from accretion discs. In section [2] we describe our 
accretion disc model, in section [3] we derive an analytical 
description for the mass supply rate from accretion discs to 
black holes which we verify in section [4] by comparison with 
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numerical calculations. In section [5] we use this analytical 
description to calculate mass outflows and apply it, as an 
example, to the mass outflow of Mrk 231 in section [S] 



2 DISC MODEL 

For our accretion disc model we assume that 

(i) the disc is rotationally symmetric, 

(ii) the disc is geometrically thin, 

(iii) the time derivative of the angular frequency is negli- 
gible, and 

(iv) the radial transport is dominated by gravitational 
and centrifugal forces. 



According to iPringlel l|l98ll ) the first three assumptions 
lead to one equation that describes the time-dependent evo- 
lution of the accretion disc's surface density E. 
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Here uj is the angular frequency, v the viscosity of the 
gas and s the distance to the central black hole, respectively. 

Using the fourth assumption and calculating the gravi- 
tational force at radius s with the monopole approximation, 
where just the disc mass within the radius s and the black 
hole mass is considered and conc entrated at s = (see, 
e.g., iMineshige fc Umemura| [T997). allows us to determine 
the angular frequency 



G (M bh + M encl («)) 



(3) 
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where 



(s) = 2tt 
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(4) 



is the disc mass that is enclosed between the inner radius 
Si n of the accretion disc and the radius s. 

As molecular viscosity is known t o be too low to ex - 
plain the observed quasar luminosities (|Lvnden- Bcll 1969), 
turbulent viscosity is held responsible for powering quasars. 
Due to the absence of a physical theory of turbulence sev- 
eral parametrizations have been deri ved to describ e the 
turbulent viscosity in accr etion discs. IShakural (|l973h and 
IShakura fc Sunvaevl l| 19731 ) parametrize it in terms of what 
became known as the a-viscosity: the product of the disc's 
scale height h, the sound speed c s and a dimensionless pa- 
rameter a which does not exceed unity. While this ansatz 
is very successful in describing non self-gravitating discs, it 
runs into serious problem s when self-gravity becomes rele- 
vant (|Duschl et aLl feoOO). We follow them and use the so- 
called /3-viscosity 
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which contains the a-viscosity in the limiting case for non 
self-gravitating dissipation-limited discs. The parameter /3 
is the inverse (critical) Reynolds number for the onset of 
turbulence of order 10~ 3 to 10 -2 . 

Solving eq. fallows us to calculate the mass supply rate 
Md from the discQ- Due to the Eddington limit, however, 
the black hole cannot necessarily swallow all the mass that 
is delivered by the accretion disc. Thus the accretion rate of 
the black hole, Mbh, is limited to 



Afhh = min M- 



(M E dd,M d ) 
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The excess mass we assume to escape from the AGN region 
resulting in the mass outflow we are interested in. 



disc's initial outer radius s out ,o leads 



M enc l,o(s) = M d , 
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As we neglect here the finite value of the inner radius Si n , this 
description is applicable only for s > Si„. But as the inner 
radius of accretion discs around black holes is of order of the 
Schwarzschild radius r s m 10~ 13 |jkk p C anc i the extension of 

accretion discs that concern us here is of order 10 2 '" 3 pc, this 
approximation does not affect our results in any significant 
way. 

We furthermore assume that the time t the gas of the 
accretion disc needs to reach its centre equals the viscous 
time-scale 



t 
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with the dimensionless parameter / of order unity. 

In the following we will discern two cases for the ratio 
between the initial enclosed disc mass and the initial central 



• In the low disc mass (LDM) case, the enclosed disc mass 
is small enough in comparison to the central object to be 
neglected for analytical approximations. This case applies 

(i) if the total disc mass is small compared to the initial 
central mass right from the beginning of the simulation, or 

(ii) for sufficiently small radii, even if the enclosed disc mass 
for larger radii is non-negligible. 

• In the high disc mass (HDM) case, the enclosed mass is 
so high that, in comparison to the central object, the latter 
may be neglected for analytical approximations. 

It is obvious that for any disc the innermost radial regions 
always belong the LDM case, while at larger radii it depends 
on the mass distribution which case applies. 



3 MASS SUPPLY RATES OF ACCRETION 
DISCS 

In the following we will derive analytical expressions for the 
mass supply rate of accretion discs. We assume an initial 
surface density profile 

So ~ s q (7) 

with the exponenlQ q g (— 2, +1) and find the disc mass 
enclosed within the radius s according to eq. [4] 
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Here we set Sj n = for simplicity. Supposing that the disc's 
initial total mass q is distributed between s = and the 



1 The square brackets in eq. [2]times (— 2ir) give the mass flow in 
the accretion disc, this e valuated at s = s; n gives the mass supply 
rate of the disc. See, e.g., IPringld jl98ll) for a detailed derivation. 

2 For q ^ —2 the integral in eq. [8] diverges, for q 1 the vis- 
cous time-scale eq. I15l is constant resp. decreasing with increasing 
radius s, which is unphysical. 



3.1 The low disc mass case 

In the LDM case the angular frequency can be derived from 
eq. [3] by neglecting the disc mass and taking into account 
the initial black hole mass Mbh.o only: 



t = 4 



with 



un = 



GMbh.o 



(11) 



(12) 



Rearranging and insertion of eq.[TT]into eq.|9] gives the total 
mass arriving at the accretion disc's centre in dependence of 
the time 



M dI (i) = M dl0 (^) 



|(9+2) 



(13) 



whose derivative is the accretion disc's mass supply rate. 

M d i(t) = M d>0 - (q + 2) f ^ij t— (14) 
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3.2 The high disc mass case 

In the HDM case the black hole mass in eq. [3] is neglected. 
With eq.[9]for the enclosed disc mass, the viscous time-scale 
is 

4(1-!) 
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We derive the mass supply rate of the disc in the same man- 
ner than above: 



, , . 2(3 + 2) 

Mdn(t) = M d ,o-^ V 1 * 1_ * • 
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3.3 The LDM/HDM transition 

Accretion from a massive disc on to the AGN's centre oc- 
curs in two phases: During an initial phase I the accreted 
material comes only from the innermost radii and thus the 
viscous time-scale for the inward movement is dominated by 
the initial black hole mass, in which case the LDM approx- 
imation applies. In the subsequent phase II, however, it is 
dominated by the disc mass, including the mass that has 
been accreted already. This is approximated by the HDM 
case. 

Equating eqs.[T4]and[l7]yields a transition time between 
the two phases: 
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/tr is of order unity as long as q is not too small (see Fig. [T]) . 



4 VERIFYING THE MASS SUPPLY RATE 
FORMULAS 

Eq. [2] describes the time- and space-dependent evolution of 
the accretion disc's surface density E. The solution of the 
set of eqs.[T]- [5] allows us to determine the mass supply rate 
Md from the accretion disc, which we will compare with our 
above analytical approach. The initial condition for the sur- 
face density E is given by eq. leaving us with the param- 
eters P (viscosity), q (initial surface density distribution), 
Sout.o (initial disc size), Md,o (initial disc mass), M^h,o (ini- 
tial black hole mass) and r\ (accretion efficiency). We set 
r\ — 10 _1 throughout this paper leaving five parameters. 

We carry out the time integration by appl ying an im- 
plicit Crank-Nicolson finite differences scheme (|Crank et al.l 
1947). The radial calculation domain extends from an in- 
ner radius Si n , which we set to 10 -4 pc for all calculations 
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Figure 1. The factors /i, fu and ftr as functions of q. 



Table 1. Parameters of the reference model. 
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presented in this paper. Tests showed that moderate varia- 
tions of this value h ave almost no effect on the res ults. Fur- 
thermore, we follow iDuschl &: Strittmatterl (|201ll l and set 
the outer radius of the computational domain to So Ut /si n . 
This is much larger than the initial physical outer radius of 
the accretion disc and minimizes the influence of the outer 
boundary condition on the temporal evolution of the model 
calculations. For numerical reasons, but with negligible in- 
fluence on the numerical results, we initially distribute 10 -3 
of the initial disc mass between the disc's outer radius and 
the numerical outer radius. Between the inner and the outer 
numerical radius 150 grid points are distributed logarithmi- 
cally. Boundary conditions are applied by setting the surface 
density to zero at the inner numerical radius which is equiv- 
alent to a vanishing torque, and by setting the mass flow to 
zero at the outer numerical radius. 

We define a reference model with the parameters given 
in Table [1] in Figs. [2](a)-(e) we vary one of the parameters 
P, q, s ou t,o, Md,o and Aibh.o from the ones of the reference 
model, respectively. We found the mass supply rate formulas 
1 141 and [TTl fit well with the numerical results for 



and 



fx 



fu 



"(9+1) 



(21) 



(22) 



with cii = 0.54 and an = 2.1. As these factors are of order 
unity (see Fig. [TJ, however, they have no crucial effect on 
our results. 

The deviations for early times in Figs. [2] (up to 10 2 yr) 
are due to setting the accretion disc's inner radius to zero 
in the analytical description. For large times (after 10 s yr) 
deviations are expected because, in contrast to the numerical 
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Figure 2. Mass supply rates of accretion discs and Eddington limit as functions of time. The dotted lines are the mass supply rates 
of the numerical calculations, the solid lines correspond to eq. 1141 and the dashed lines to eq. 1171 The crossed line is the Eddington 
limit (Fig. (f) only). In Fig. (a) we varied the /3-parameter, in Fig. (b) the exponent q, in Fig. (c) the disc's initial outer radius s O ut,0) in 
Fig. (d) the initial disc mass M<j q and in Fig. (e) the initial black hole mass Af^h o from the reference model, respectively. Fig. (f) shows 
numerical and analytical mass supply rates and the Eddington limit for our model parameters of Mrk 231. 
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results, eq. [9] does not take into account that the disc has a 
finite radius at s ou t o- 



5 MASS OUTFLOWS 

For the following order-of-magnitude calculations we use an 
average mass supply rate for phase II, rather than the one 
derived in eq. [17] 



,Vavg _ M d ii(t) 
MdU ~ t 



2(g + 2) 

= M d , (—) t fri 



(23) 



The average rate for phase I may be calculated in a similar 
manner. 

As we will show in section [6] for not too large times and 
for disc masses much larger than initial black hole masses the 
mass supply rate of the disc by far exceeds the Eddington 
limit 



M d > M Edd • 



(24) 



so almost all mass delivered by the disc escapes from the 
AGN and the excess mass flow, i.e., the difference between 
what is delivered by the disc and the Eddington limit, can 
be approximated as 



M d ,0 
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(25) 



This mass we assume to escape from the system as an out- 
flow, the lifetime of which can be estimated as 



tit 



R 



(26) 



with R the extension and v the velocity of the outflow, re- 
spectively. Thus eq. [25] gives 



3(q+l) 

•<77) ~ q U -l./ J 



2(g + 2) 
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(27) 



which connects the parameters of the outflow, Af cxc , R and 
v, with the parameters of the accretion disc, /3, q, s ou t,o and 
M d , ._ 

Af e xc can be determined by measuring the velocity, the 
extension, the solid angle and the column density of the out- 
flow. For different approaches to determine M e xc fro m these 
quantities, see, for instance, IShih fc Rupkel (|201Ch (thin- 
shell model) and iRupke et al.l l|2005l ) (thin- and thick-shell 
model) . 



of the outflow of order ~ 3kpc, a velocity of ~ 10 kms -1 
and a mass outflow rat^fl 420 Mq yr -1 . 
With these data we write eq. [23 as 

1-9 



M c: 
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with 



g(q) = 9.85 • 2.43" ■ 10 



-10-79 



fix 
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(28) 



(29) 



For Mrk 231 the brackets in eq. [28] equal unity and 
by solving the equation g(q) = 1 we can determine q « 
— 1.38 and thus the (time-averaged) surface density profile 
in Mrk 231 as E oc s" 138 . 

The corresponding mass flows feqs. I14l and ll7[) we plot 
in Fig. [2](f) and compare them with numerical solutions of 
this model to show the agreement of numerical and ana- 
lytical description also for these parameters. Additionally 
Fig. [21(f) shows that the Eddington limit is indeed much 
smaller than the disc's mass supply rate for not too large 
times as we assumed in the previous section. 

The transition time tu ~ 300 yr is small compared to 
the lifetime of the outflow of tu » 3 • 10 6 yr, which justifies 
the use of the phase II model for describing the mass outflow 
from Mrk 2310 



7 SUMMARY 

We derived analytical expressions for the mass supply rate 
of accretion discs which depend on basic disc parameters, 
black hole mass and time and verified them by comparison 
with the results of numerical calculations. 

When the Eddington limit is much lower than the mass 
supply rate of the accretion disc, this allows us to constrain 
the properties of mass outflows. 

With these analytical expressions we are able to connect 
the basic properties of the outflow (M oxc , R and v) and 
of the accretion disc (/?, q, s ou t,o, Md,o)- This relationship 
can be used to estimate basic properties of quasars. As an 
example, we applied it to Mrk 231 and found the surface 
density distribution of its accretion disc to be oc s -1,38 . 

As in our scenario the outflow is driven by the limitation 
of the black hole's accretion rate through the Eddington 
limit, we conclude that, at least in such cases, (radiation) 
pressure is capable of driving huge quasar mass outflows. 



6 MASS OUTFLOW FROM MRK 231 

Mrk 231 is classifi ed as ultra-luminous infrared galax y 
IIRieke fc Lowll 19721 ) and Seyfert I galaxy |Weedmanll 19731) . 
It hosts an accretion disc with a mass of 3.15 ■ 10 10 Mq , 
an outer radius of 1.7 kpc l|Downes fc Solomon! 1 1998 | ) and 
a central unresolved mass of < 8 • 10 6 Mm l|Richards et al.1 
2005), which we assume to be the central black hole's mass. 
For the /3-parameter we adopt the value /3 = 0.005 from our 

refe rence model. 

IRupke fc Veilleuxl l|201ll ) recently described a massive, 
wide-angle outflow from Mrk 231 and estimate an extension 
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